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ABSTRACT 

We present deep Hubble Space Telescope Wide Field Planetary Camera 2 (WFPC2) optical ob- 
servations obtained as part of the ACS Nearby Galaxy Survey Treasury (ANGST) as well as early 
release Wide Field Camera 3 (WFC3) ultra-violet and infrared observations of the nearby dwarf star- 
bursting galaxy NGC 4214. Our data provide a detailed example of how covering such a broad range 
in wavelength provides a powerful tool for constraining the physical properties of stellar populations. 
The deepest data reach the ancient red clump at VIfsiaw ~-0.2. All of the optical data reach the 
main sequence turnoff for stars younger than ^300 Myr, and the blue He burning sequence for stars 
younger than 500 Myr. The full CMD-fitting analysis shows that all three fields in our data set are 
consistent with ^75% of the stellar mass being older than 8 Gyr, in spite of showing a wide range in 
star formation rates at the present day. Thus, our results suggest that the scale length of NGC 4214 
has remained relatively constant for many Gyr. As previously noted by others, we also find the galaxy 
has recently ramped up production, consistent with its bright UV luminosity and its population of 
UV-bright massive stars. In the central field we find UV point sources with F336W magnitudes as 
bright as -9.9. These are as bright as stars with masses of at least 52-56 M Q and ages near 4 Myr 
in stellar evolution models. Assuming a standard IMF, our CMD is well-fitted by an increase in star 
formation rate beginning 100 Myr ago. The stellar populations of this late-type dwarf are compared 
with those of NGC 404, an early-type dwarf that is also the most massive galaxy in its local environ- 
ment. The late-type dwarf appears to have a similar high fraction of ancient stars, suggesting that 
these dominant galaxies may form at early epochs even if they have low total mass and very different 
present-day morphologies. 

Subject headings: galaxies: individual (NGC-4214) — galaxies: stellar populations — galaxies: spiral 
— galaxies: evolution 



1. INTRODUCTION 

Irregular galaxies are the dominant type of star form- 
ing galaxy by number. Despite their typically isolated 
environments and relatively low mass, the star forma- 
tion histories (SFHs) of irregular galaxies are surpris- 
ingly complex and varied. In the context of hierarchical 
structure formation, the formation epoch of these galax- 
ies is expected to be younger than that of more massive 
galaxies, a phenomenon referred to as cosmic downsizing. 
Deep resolved stellar photometry allows us to measure 
their formation epoch, and to examine their recent SFHs 
in detail. 

Historically, early-type dwarfs were thought to be 
dominated by an cient populations ( Baad e~fc Gaposchkinl 
119631 : fZinn 1980), while dlrrs have been considered to be 
the result of continuous, relat i vely constant , star forma- 
tion (e.g. IHunter et all 119821: IHunterl Il997l ). This pic- 
ture has broken down in recent decades. Early-type 
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dwarfs were found to har bor significa nt populations of 
intermediate-age stars (see|Mateo 1998, for a review, and 
many other studies), and recent studies of dlrrs strongly 
suggest that they u ndergo episodic star formation (e.g. 
Stin son et al.ll2007l and references therein). 

If dlrrs are characterized by relatively constant star 
formation, we might expect the fraction of ancient stars 
to be significantly lower in dlrrs compared to early- 
type dwarfs. However, resolved stellar photometry of 
local dwarfs has shown that the fraction of ancient 
stars is roughly the same in both m orphological types 
(jWeisz. D. and Williams et al.l l2010f ). Other resolved 
stellar population studies suggest bursts of about 500 
Myr dura tion play an importan t role in the evolution 
of dwarfs ()McQuinn et al.l [2010al fbh . It is therefore pos- 
sible that the presence of recent star formation episodes 
may be the only difference between these classes of dwarf 
galaxies and that such recent episodes of star formation 
actually contribute very little to the total stellar mass 
of the systems. Therefore, the difference between these 
distinct morphological types may be due to only a very 
small fraction of the stellar mass. To further investigate 
the possibility that the vast majority of the stellar mass 
in dlrrs can be as old as the stellar mass of early type 
dwarfs and that the star formation events witnessed in 
dlrrs are episodic, we focus our attention on the nearby 
galaxy NGC 4214. 

NGC 4214 is classified as an irregular galaxy (IAB(s)m, 
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Ide Vaucouleurs et al.ll 19911) with a mass of ^1.5xl0 9 M Q 
(jKarachentsev et al J 120041) . similar to the Large Magel- 
lanic Cloud. It is one of the nearest example s of a star- 
burst , and was the hos t of a Type I supernova (iWellmannl 
119551) and recent nova ([Humphreys et al. 2010). It there- 
fore has been studied in detail, especially at UV and NIR 
wavelengths that probe the you ng stellar populations. 
For example, iFanelli et al.l (|1997|) were able to make use 
of far-UV imaging from the Ultraviolet Imaging Tele- 
scope and ground-based I-band to decompose the galaxy 
into a centrally concentrated intense star forming com- 
ponent and an extended smooth disk. They suggest the 
two component structure is the result of a recent merger. 

The bright resolved stellar populations have also been 
characterized. A dominant red giant branch was de- 
tected, and a larg e meta llicity spread was measured by 
IDrozdovskv eTalT (I2002T) using WFPC2 and NICMOS 

photometry. lUbeda et al.l (l2007allbT l used HST/WFPC2 
and HST/STIS imaging to suggest that the stellar ini- 
tial mas s funct ion appears to be steeper than that of 
iSalpeterl (|1955l ) at high masses (>20 M ), even though 
the galaxy's metallicity is similar to that of the Small 
Magellanic Cloud. However, it is unclear how this slope 
is affected by fluctuations in the star formation rate. 

We further examine the stellar populations of 
NGC 4214 using deeper data than previous studies. The 
new data include two HST WFPC2 fields of the outer 
disk and one WFC3 field of the inner regions. Through 
detailed stellar evolution model fitting, we measure the 
star formation history (SFH) of these regions and look for 
differences between quiescent and intensely star forming 
regions. Section 2 describes our data set and photometry. 
Section 3 presents the results of our model fitting proce- 
dure. Section 4 discusses the results of the measurements 
in the context of galaxy morphology and environment, 
and Section 5 summa rizes our conclusions. We assume 
(m-M)o = 3.03 Mpc (|Dalcanton et al.ll2009D for conver- 
sions of angular measurement s to physical distances and 
assume an inclination of 38° (Karachen tsev et al.l 120041) 
for surfa ce density measurem ents. We adopt a five-year 
WMAP (iDunklev et al.1 120091 ) cosmology for all conver- 
sions between time and redshift. 

2. DATA 

As part of the ANGST program after the failure of 
ACS, from 2007-Dec-04 to 2007-Dec-23, we performed 
deep WFPC2 observations a field in the NGC 4214 
disk located at R.A. (2000)=183. 847375 (12:15:23.4), 
deck (2000)=36.362333 (+36:21:44) with a rotation an- 
gle PA_V3=119.9 degrees (GO-10915). To improve our 
radial coverage to be closer to the goal of the origi- 
nal ANGST program, from 2009-Feb-23 to 2009-Feb-26, 
we performed shallower supplemental observations for 
a field located at R.A. (2000)=183.878103 (12:15:30.7), 
deck (2000)^36.356808 (+36:21:24.5) with a rotation 
angle PA_V3=51.0 (GO-11986). At the same time, 
WFC3/IR and WFC3/UVIS imaging of the central area 
of the galaxy, proposed by the Scientific Oversight Com- 
mittee (SOC), were released. Both fields were observed 
from 2009-Dec-22 to 2009-Dec-23 located at R.A. (2000) 
=183.913333 (12:15:39.2), deck (2000) = 36.326944 
(+36:19:37.0) with a rotation angle PA_V3=120.0 (GO- 
11360, PI: O'Connell). Figure [T] shows outlines of the 
fields' locations. Our field locations were chosen to max- 



imize the number of disk stars and avoid crowding. 

In the deep field, we obtained 15 full-orbit exposures 
with the WFPC2 through the F606W (wide V) filter, 
and 29 full-orbit exposures through the F814W (J equiv- 
alent) filter. These data totaled 39000 s and 75400 s of 
exposure time in F606W and F814W, respectively. In 
the other WFPC2 field, we obtained 2 orbits through 
F606W, totaling 4800 s, and 4 orbits through F814W, 
totaling 9600 s. The WFC3/IR data contained 1198 s 
and 2398 s of exposure in F110W and F160W, respec- 
tively. The WFC3/UVIS data contained 1683 s, 1540 s, 
and 1339 s in F336W, F438W, and F814W, respectively. 
All WFPC2 images were calibrated in the HST pipeline 
with CALWP2 using OPUS version 2006_6a for the 2007 
data and 2008_5c for the 2009 data. All WFC3 images 
were calibrated in the HST pipeline with CALWF3 ver- 
sion 2.0. 

All of the data were processed through the ACS 
Nearby Gal axy Survey Treasury (A NGST) data analy- 
sis pipeline (|Dalcanton et akll2009D. updated to i nclude 
WFC3 UVIS and IR data (IDalcanton et akl[20Tl . As a 
brief summary, the photometry was measured simulta- 
neously for all of the objects in the uncombined images 
using the software p ackages HSTPHOT (for WFPC2) 
and DOLPHOT 1.2 (|Dolphinll200"ol 120091) including the 
WFC3 module. These packages are optimized for mea- 
suring photometry of stars on HST images using the well- 
characterized and stable point spread function (PSF) 
calculated with TinyTim0 The software fits the PSF 
to all of the stars in each individual frame to find PSF 
magnitudes. It then determines and applies the aperture 
correction for each image using the most isolated stars, 
corrects for the charge transfer efficiency of the WFPC2 
detectoiQ, combines the results from the individual ex- 
posures, and converts the measured count rates to the 
VEGAmag system. 

Our photometric error and completeness were then as- 
sessed by running at least 10 6 artificial star tests for each 
field, in which a star of known color and magnitude was 
placed into the images and the photometry rerun to de- 
termine if the star was recovered and, if so, the difference 
between the input and output magnitude. 

The photometry output was then filtered to only 
allow objects classified as stars with signal-to-noise 
>4 in both filters. The list was further culled us- 
ing sharpness {\F606W sharp + F8UW sharp \ < 0.27 
for WFPC2, \F336W sharp + F438W sharp \ < 0.27 for 
UVIS, and \F110W sharp + F160W shorp | < 0.35 for IR) 
and crowding (F606W crowd + F%lAW crowd < 0.7 for 
WFPC2, F336W crowd + F438W crowd < 0.7 for UVIS, 
and F110W crowd + Fl60W crowd < 0.48 for IR). The 
sharpness parameter returns zero if a star is perfectly 
fit, negative if it is b roader, and po sitive if it is sharper 
than a typical star (|Dolphinl I2000D . The crowding pa- 
rameter gives the difference between the magnitude of 
a star measured before and after subtracting the neigh- 
boring stars in the image. When this value is large, it 
suggests that the star's photometry was significantly af- 
fected by crowding, and we therefore exclude it from our 
catalog. The values for the sharpness and crowding cuts 

5 http:/ /www. stsci.edu/software/tinytim/ 

6 July, 2008 formulae 

http: / /purcell. as. arizona.edu/wfpc2_calib/2008_07_19. html 
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for WFPC2 were the standard used by the ANGST pro- 
gram. The cuts for WFC3 were chosen by looking at 
CMDs resulting from different possibilities and choos- 
ing values that produced clean CMD features and a low 
number of contaminants, falling outside of any known 
feature. We note that our tests of different cuts showed 
that the choices of cuts have little impact on our final 
SFH measurements outside of the most recent time bin 
(<10 Myr), most likely due to the significant clustering 
of the very youngest stars. Our final optical star catalogs 
contained 16806, 28088, 138350, and stars for the outer 
(WFPC2 N), bridge (WFPC2 NW), and central (UVIS) 
fields, respectively. 

Our final UV catalog for the central field contained 
77656 stars, and our final IR catalog for the central 
field contained 19807 stars. In the end, our analysis of 
the optical (F438W and F814W) photometry from UVIS 
yielded results that were consistent with those of the UV 
(F336W and F438W) photometry from UVIS. Since no 
additional insight was gained from the optical photom- 
etry, we will not discuss it further but include the SFH 
and limiting magnitudes for reference. Any statements 
made about the UV photometry from UVIS was con- 
firmed with the optical photometry from UVIS. The final 
CMDs are shown in Figure [5J 

2.1. Field Division 

While all the WFPC2 imaging of the outer portions of 
the disk provided high-quality photometry going much 
fainter than the tip of the red giant branch (and even 
fainter than the red clump in the deep field), the central 
WFC3 fields were not as consistent. The UVIS images 
proved to have similar completeness limits throughout 
the field, but the IR images were clearly over-crowded 
near the galaxy center. This strong differential crowding 
in the IR data made analysis of the full IR field dubi- 
ous. We therefore separated the innermost portion of 
the central field into two spatial components, following 
the isophotes of the galaxy. Our division ellipse and the 
galaxy isophotes in the F110W image are shown in Fig- 
ure 131 

3. CMD FITTING AND TEMPORAL RESOLUTION 

We measured the star formation rate and metallicity 
as a func tion of stellar age using the software package 
MATCH ()Dolphinll200"l to fit the observed CMDs. We 
adopted magnitude cuts set to the limits provided in 
Table [H and t hen fi tted the stellar e volution models of 
iGirardi et all (120021 wi th updates in lMarigo et all 120081 
and lGirardi et al.1 12008). convolved with our photometric 
error and comp leteness statistics and populated with a 
ISafpeterl ([1955) initial mass function (IMF). The choices 
of software and model s used for the ANGST p roject are 
discussed in detail in Williams et al.l ([2009b). Briefly, 
comparisons of results when fitting CMDs using differ- 
ent models and fitting software have shown that the re- 
sults are consistent within the estimated uncertainties; 
furthermore, our choice of models provides the largest 
range of ages and metallicities publicly available. 

We first fit the data assuming a single foreground red- 
dening Ay=0.07 and distance m — Mp=27.41, which 
were obtained from the ANGST sur vey ([Dalcanton et al.l 
2009), and the lSchlegel et al.l (jl998l ) Galactic dust maps. 
The best fit provides the relative contribution of stars of 



each age and metallicity in each field. For the shallowest 
data, the IR and optical data for the innermost region, 
we limited the number of free parameters by imposing 
an "increasing metallicity" constraint on the fit, whereby 
the metallicity of the population was not allowed to de- 
crease with time (within the measured errors) . These fits 
are turned into SFHs of star formation rate and metal- 
licity as a function of time. 

To assess the uncertainties of our fits, we then ran 
Monte Carlo fitting tests. These tests assess 2 types 
of errors: random errors due to Poisson sampling and 
errors in photometry, and systematic errors due to de- 
ficiencies in the stellar evolution models as well as any 
offset in distance, reddening, and/or zero-points. The 
Poisson errors are accounted for by resampling the con- 
volved best-fit model 100 times. Then, when fitting each 
of these realizations, the systematic errors are accounted 
for by introducing small random shifts in the bolomet- 
ric magnitudes and effective temperatures of the models. 
These shifts are introduced at the level of the differences 
between models in the literature, and therefore serve as 
a proxy of the effects of our choice of stellar evolution 
models. 

Our final uncertainties were calculated as the 68% con- 
fidence intervals of the results of all of our Monte Carlo 
test fits. These total uncertainties are used as the error 
bars in all subsequent plots and analysis. 

We can also assess the degree of systematic errors us- 
ing plots of residuals between the observed CMD and 
the model fit. An example of the residuals of our fits is 
shown in Figure 2J Our fits assume the Galactic fore- 
ground extinction over the field. Therefore, the low-level 
residuals over the entire CMD for these fits suggests that 
our photometry is dominated by stars in front of any sig- 
nificant internal dust layer in the NGC 4214 disk, with 
a modest amount of internal extinction (cf. § 14. 2p caus- 
ing the slight excess in the data redward of the model 
main-sequence. Furthermore, the lack of any strong fea- 
tures in the fit residuals suggests that the distribution 
of stellar luminosities is well-fitted assuming a Salpeter 
IMF, once the effects of the age distribution are modeled. 
These fits show that the main sequence luminosity profile 
interpreted bv lUbeda et al.l (2007a. b) as suggesting vari- 
ations in the IMF, is also consistent with a standard IMF 
with a varying star formation history; thus IMF varia- 
tions and SFH are somewhat degenerate in such a com- 
plex region. Our uncertainties come from measurements 
assuming a constant IMF. Therefore, our uncertainties 
would increase if variations in the IMF were allowed. 

We used a very fine grid of model stellar isochrones 
(0.05 dex in age and 0.1 dex in [Fe/H]) to provide the 
best possible fit to our data. Because the grid is very 
fine, pushing star formation between adjacent time bins 
has little effect on the quality of the resulting fit. The 
degree to which this degeneracy is true depends on the 
quality of the data being fitted, which is characterized 
by our Monte Carlo tests. 

The uncertainties in the cumulative age distribution 
cannot be improved by further binning of the data in 
time. Therefore, we plot the full time resolution of our 
fit to the cumulative age distribution, showing the un- 
certainties from our Monte Carlo tests, which include 
the effects of pushing star formation between adjacent 
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time bins. 

On the other hand, within a given time bin, the un- 
certainty on the measured mean SFR is related to the 
sensitivity of the CMD to stars in that age bin. As the 
size of the bin increases, more locations of the CMD will 
be affected by the bin. Thus, increasing the bin size in- 
creases the sensitivity of the CMD fit to that time bin, 
thereby reducing the uncertainty on the mean rate within 
that time bin. It is therefore possible to reduce the un- 
certainties on the measured SFR by increasing the length 
of the time bins. 

To determine optimal time bins for our differential 
SFH, we apply the uncertainties in our cumulative age 
distribution, as determined from our Monte Carlo tests. 
Specifically, we define bins for which the cumulative 
fraction of stars increased with statistical significance. 
Specifically, the bins are defined so that at each bin 
boundary, the la upper-limit of the cumulative fraction 
of stars formed in the older bin must be less than the la 
lower-limit of the younger bin. Thus, each bin contains 
enough signal in the CMD that if the bin were removed, 
the fit to the data would be significantly degraded. We 
note that for the optical and UV fits, we adopted the 
time bins measured from the IR fits because the optical 
and UV data were not sensitive enough to the old stellar 
populations to provide cumulative distributions. 

4. RESULTS 
4.1. Star Formation History 

The resulting SFHs for our CMD fits to the 4 regions 
are shown in Figure [5] The UVIS data did not pro- 
vide sufficient depth to probe ages greater than ^200 
Myr, and thus no points are plotted for older ages. Cu- 
mulative plots of the stellar mass formed are shown in 
Figure [51 and the results for all regions are overplotted 
in Figure [7J Since the UVIS data were not sensitive to 
old populations, we could not generate cumulative SFHs 
from the UVIS data. 

The SFHs show that in all of the areas observed, the 
majority of the stellar population is old. In the cen- 
tral portion of the galaxy, the relatively shallow data 
are only able to constrain the population only weakly, 
putting a lower limit of >80% of the stellar mass having 
ages >4 Gyr. In the outermost (and deepest) field, the 
constraint is tighter, with >74% of the stellar mass hav- 
ing ages >8 Gyr. Taking uncertainties into account, all 
of our data are consistent with this constraint from the 
deepest photometry. 

The consistency of all of our data, from the optical to 
the IR, from the center to the outskirts of the disk, with 
more than three quarters of the stellar mass being very 
old is surprising. This galaxy is well-known for being a 
young s tarburst galaxy, being, know n as a "Wolf-Rayet" 
galaxy ([Sargent fc Filippenkolll99ll) . While the current 
star burst is already known to be responsible for onl y 
a few percent of the stellar mass (|Huchra et al.lli~983l ). 
it is still surprising that such a small percentage of the 
population is young (<1% in the past 50 Myr, <4% in 
the past Gyr), and such a high percentage is very old. 
Overall, the galaxy has evolved very little since z ~ 1, 
as shown in our plot of the surface density profile over 
the past 6.3 Gyr (Figure [8]). 

Our measurements also provide some indication of 



the metallicity of the stellar populations in NGC 4214. 
Overall, the metallicity of the populations fell in the 
range of -1.6<[M/H]<0, with the old population being 
-1.6<[M/H]<-0.6 in the IR data and -1.4<[M/H]<-0.6 
in the WFPC2 data, and the young population being - 
0.5<[M/H]<0 in the UVIS data. The WFPC2 data had 
very little sensitivity to metallicity for the young pop- 
ulation, as MS stars of different metallicities have very 
similar colors in F606W-F814W. 

4.2. The Young Stellar Population 

The old population likely dominates at all radii, sug- 
gesting that this population is well- mixed and NGC 4214 
is old. However, the young stellar populations are clearly 
not well-mixed, as their density clearly varies signifi- 
cantly with position in the galaxy. This difference can be 
seen by looking at the strong core He-burning sequence 
extending vertically from the tip of the RGB in the cen- 
tral WFC3/IR CMD and northern WFPC2 CMD but 
relatively absent in the northwest WFPC2 CMD. The 
central starburst is evident, even when the population 
is probed only in the IR, where massive young stars are 
faint, although we note that the very youngest popula- 
tion (<10 Myr) in the innermost region is not reliably 
measured in the IR (Figure [5]) . The SFH of the inner- 
most region of the galaxy is the only one with a sustained 
recent ( ;$ 30 Myr) star formation rate that is higher than 
any rate in previous epochs. This result does not mean 
that the rate is higher than it has ever been in the past, 
since our time resolution would not allow us to detect 
short, strong bursts at ancient times, but the result does 
show that very recent star formation near the galaxy cen- 
ter is at least a factor of 2 higher than the average rate 
over the galaxy's history. Just outside of this central 
starburst, the star formation rate is slightly less than the 
overall mean rate but higher than the mean rate over the 
past several Gyr. 

Outside of the central WFC3 field, the knots of star for- 
mation in the northern WFPC2 field are easily seen in the 
CMD. Recent star formation produces both the strong 
red He-burning sequence extending brightward from the 
RGB and the well-populated upper MS. The SFH shows 
a strong increase in star formation rate starting 100 Myr 
ago with a prominent peak at very recent times (5-12 
Myr ago), providing some indication of the age of these 
features. Prior to 100 Myr ago, our analysis shows that 
these regions of NGC 4214 were relatively quiescent, per- 
haps bearing resemblance to an early-type dwarf. 

In Figure |H] we plot the spatial distribution of young 
stars as a function of brightness. The brightest 
(youngest) stars appear to be more tightly clustered than 
the fainter stars. A 2-dimensional K-S test yields only 
a 4% probability that the brightest MS stars (cyan and 
white diamonds in Figure EJ and the faint MS stars (red 
diamonds in Figure are drawn from the same spa- 
tial distribution. Thus stars have likely migrated signif- 
icantly from their birth locations on timescales of just 
~12 Myr (the main-sequence lifetime of a 12 M Q star). 

To gain a sense of the velocities with which stars are 
spreading from their birth regions, we looked for a good 
example region to study in detail. One of the knots of 
star formation (RA = 183.8884, DEC = 36.3683, see Fig- 
ure [TO]) was isolated enough to look for a correlation be- 
tween the magnitude of the brightest MS stars and the 
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distance from the center of the knot. Taking the bright- 
est MS star in annuli of 0-6", 6-8", 8-10", and 10-12" 
(see Figure [TO]) , there is a weak correlation between the 
magnitude of the brightest MS star and distance from the 
central knot. Each 75 pc farther out, the brightest MS 
star is — 1 magnitude fainter. If we assume the bright- 
est MS star is a proxy for the MS turnoff magnitude, 
at these bright magnitudes on the upper main sequence 
1 magnitude corresponds roughly to 10 Myr of age. If 
we further assume that the weak correlation is real and 
that it is due to the diffusion of stars produced in the 
center of the knot, the slope of the correlation is con- 
sistent with a diffusion speed of —8 pc Myr -1 (—8 km 
s _1 ). This value is similar to t he velocity dispersion o f 
B stars in the Milky Way disk (jDehnen fc Binnevi[T998l) 
and t he diffusion speed foun d for an outer spiral arm in 
M81 (|Williams et al.ll2009bl) . 

The brightest MS stars in our catalog provide diagnos- 
tic tests for massive stellar evolution models. HST UV 
spectroscopy has been measu red for the central star burst, 
yielding an age of 4-5 Myr (jLeitherer et al.lll996D . Our 
photometry reveals point sources with F336W ([/-band 
equivalent) magnitudes of —17.6, which is Mp336W=- 
9.9 at the distance and extinction of NGC 4214. We 
show our F336W-F438W CMD in Figure CU along with 
the isochrone for an age of 4 Myr and a metallicity of 
[M/H]=-0.4 shifted to the distance and foreground ex- 
tinction of NGC 4214. Assuming these objects are not 
unresolved compact clusters, these bright stars are con- 
sistent with the brightest blue portion of this isochrone, 
which represents model stars with masses of 52-56 M . 
The low foreground extinction (Ay — 0.07) and good 
agreement between the data and model isochrone sug- 
gests that these stars are not strongly affected by extinc- 
tion internal to NGC 4214. NGC 4214 is known to con- 
tain just a modest amount of internal extinction (inter - 
nal GALEX FUV extinction -0.58 mag. lLee et "aT1[2009l) . 
and these particular UV-bright stars are apparently be- 
tween us and any of significant dust within NGC 4214, 
making them easily detected and measured in the UV. 

Finally, we estimated the approximate total SFR for 
the galaxy over the past —100 Myr, as it has been rela- 
tively high during most of this period. If we add together 
the SFRs from our 4 regions, we obtain a total rate of 
-0.1 M yr" 1 , -80% of which is in the central UVIS 
field. This rate is consistent with the SFR recently mea- 
sured for NGC 4214 fro m HST/WFPC2 photometry by 
IMcQuinn etaD (|20l0allbl ) . 

4.3. Dominant Galaxies in their Environment 

Interestingly, the overall population is similar to that 
of the very different nearby dwarf galaxy NGC 404, 
which is a somewhat m ore massive (— 4.5xl0 9 M ; 
iKarachentsev et al.1 120041 ) SO galaxy with very little re- 
cent star formation. The ANGST data for that galaxy 
show that outside the inner regions —90% of the stel- 
lar population is >8 Gyr old and —75% is >10 Gyr old 
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(|Williams et al.l 12010ft . While the overall stellar popu- 
lations of NGC 404 are older than those of NGC4214, 
the difference is not overwhelming. In fact, by 4 Gyr 
ago, both galaxies had formed 90% or more of their stel- 
lar mass. Apparently, the vast difference in morphology 
between these 2 galaxies is mainly due to just a small 
percentage of the stellar mass. Our results show that 
—1% of the stellar mass in NGC 4214 formed in the past 
—100 Myr. The percentage of stellar mass formed in the 
past 100 Myr is much smaller in NGC 404, only <0.1% 
outside of the galaxy center. However, there is evidence 
for a 1 Gyr old starburst in the center of NGC 404 
(|Seth et aljl2010h : about 1 Gyr ago NGC 404 may have 
more closely resembled NGC 4214's present-day appear- 
ance. While this difference represents a small percentage 
of the total stellar populations of the galaxies, it results 
in their vastly different morphological classifications. 

Both NGC 4214 and NGC 404 appear to have older 
median ages than some more massive disks, such as M33 
(| Williams et al.l l2009al) . Perhaps their old stellar ages 
are somehow due to their relationship with their envi- 
ronments. Both galaxies are the most massive in their 
local environment. Perhaps it is common for the domi- 
nant galaxy by mass within any group of galaxies (where 
NGC 404 is isolated, making its own "group") to be dom- 
inated by ancient stars. If this effect is indeed common, 
it may indicate that the most massive galaxy in a group 
dominates gas accretion in the early stages of evolution, 
winning the local "downsizing" battle to be the dominant 
star forming object during the epoch of formation. 

5. CONCLUSIONS 

We have analyzed deep HST/WFPC2 photometry of 
the NGC 4214 disk and HST/WFC3 photometry of the 
central portion of NGC 4214. Full CMD modeling of the 
photometry shows that the stellar populations through- 
out the disk are old, with —75% of the stellar mass 
older than —8 Gyr. This result shows that overall the 
stellar populations of NGC 4214 are similar to those of 
the SO galaxy NGC 404, other than the youngest 1% of 
the stellar mass, which formed in the past 100 Myr in 
NGC 4214. The similarity suggests that a few hundred 
Myr ago, NGC 4214 may have looked very similar to 
NGC 404 today, though NGC 4214 may have been more 
gas- rich. Alternatively, NGC 4214 may have recently ac- 
quired gas fr om a merging satellite, such as the merger 
suggested bv iFanelli et ail (|1997[ ). NG C 4214 currently 
has m ore than 10 nearby satellites (e.g. iDalcanton et al.1 
2009, Fi gure 2); therefore such an alternative is a strong 
possibility. Finally, we presented an argument that be- 
cause both galaxies are the dominant members of their 
local environments, they were in favorable positions to 
form such a high percentage of their stars so early. 
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TABLE 1 

Magnitude Limits of the Observed Regions 



Region 


F336iy 50 a 


F438W50 b 


F606W 50 C 


F814W 50 d 


FU0W 50 c 


F160Wso f 


GCDS 


UV INNER 


26.0 


27.1 










0.34 


UV OUTER 


26.1 


27.2 










1.15 


IR INNER 










23.5 


22.7 


0.52 


IR OUTER 










24.0 


23.2 


0.99 


BI INNER 




26.9 




26.0 






0.48 


BI OUTER 




27.1 




26.2 






1.10 


WFPC2 N 






26.7 


26.1 






2.27 


WFPC2 NW 






27.9 


27.1 






3.17 



a The 50% completeness limit of the F336W data. 
b The 50% completeness limit of the F438W data. 
c The 50% completeness limit of the F606W data. 
d The 50% completeness limit of the F814W data. 
e The 50% completeness limit of the F110W data. 
f The 50% completeness limit of the F160W data. 

g The median dcprojectcd galactocentric distance of the stars in the region in kpc. 




Fig. 1. — The locations of our NGC 4214 fields are shown on a 3-color image using Sloan Digital Sky Survey u' (blue), g' (green), and %' 
(red) images. North is up. East is left. The small central field is WFC3/IR. The larger central field is WFC3/UVIS. The outer fields are 
WFPC2, with the westernmost being the deepest. 




Fig. 2.— The CMDs of our 4 fields (§|2t. Top Left: WFC3/IR. Top Right: WFC3/UVIS. A single isochrone for an age of 4 Myr and a 
metallicity of [M/H]=-0.4 is overplotted f § 14.21 1, Arrow indicates direction and magnitude of Ejy— b=0.3. Bottom Left: the northern WFPC2 
field. Arrow indicates direction and magnitude of £f606W-F814W = 0-3- Bottom Right: the northwestern WFPC2 field with isochroncs 
overplotted (from blue to red: [M/H]=-0.4 and log(age) = 7.3,7.6,8.0,8.3,8.6, followed by log(age)=10.0 and [M/H] = -1.7,-1.3,-0.7,-0.4, 
respectively). Arrow indicates direction and magnitude of i?F606W-F8l4W=0.3. 
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Fig. 3.— A WFC3 image of the central 2'X2' of NGC 4214 shown in a 3-color image using WFC3 F438W (blue), F110W (green), and 
FWOW (red) images. North is up. East is left. Isophotal contours are shown in green. Our field division ellipse is shown in white. The 
inner region is much more crowded and dusty than the outer region, making CMD fitting more reliable when the two regions are modeled 
separately. 
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Fig. 4. — The observed (upper-left), best-fit model (upper-right), difference (lower-left), and difference normalized by uncertainty (lower- 
right) CMDs of our fit to the UVIS data of the inner region of NGC 4214. Scale bars are shown to the right of each plot. The fit is generally 
good over the entire CMD assuming foreground extinction, but also shows a modest amount of internal extinction (cf. § 14,211 which causes 
the slight excess in the data redward of the model main-sequence. 
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Fig. 5. — The SFHs of our fields. Upper left: The SFH of the WFPC2 field to the north of the center of the galaxy. Upper right: The 
SFH of the deep WFPC2 field to the northwest of the center of the galaxy. Second row left: The SFH of the inner region of the WFC3/IR 
field. Second row Right: The SFH of the WFC3/IR field excluding the inner region. Third row left: The SFH of the inner region of the 
WFC3/UVIS field from the F438W and F814W photometry. Third row right: The SFH of the WFC3/UVIS field excluding the inner region 
from the F438W and F814W photometry. The SFH for times prior to ~200 Myr ago is unconstrained by these data. Bottom left: The SFH 
of the inner region of the WFC3/UVIS field from the F336W and F438W photometry. Bottom Right: The SFH of the WFC3/UVIS field 
excluding the inner region from the F336W and F438W photometry. The SFH for times prior to ~200 Myr ago is unconstrained by the 
UVIS data. 
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Fig. 6. — The cumulative formation of stellar mass calculated from the SFHs of our fields. The central WFC3 field has been broken up 
into an inner region and an outer region, shown in the top two panels. The UVIS data did not constrain the cumulative formation because 
it contained no significant information concerning stars older than 150 Myr. 




Fig. 7. — All cumulative formation histories overplotted. Darker shades of gray represent overlaps in the cumulative SFHs at different 
radii. All of our measurements are consistent with 75% of the stellar mass forming by z~l. The best-fit solutions are shown with colored 
lines (orange=N4214_NW.WFPC2, green=N4214_N_WFPC2, blue=N4214_OUTER_IR, purple=N4214_INNERJR). 
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Fig. 8. — The cumulative stellar mass density profile of NGC 4214 at 3 epochs, as calculated from our SFH measurements. Changes in 
the IMF and/or mass cutoff could shift all of the surface density data up or down. The scale-lengths are unaffected by these IMF issues. 




Fig. 9. — The locations of the massive (>12 Mg) young stars in our images. Bluer colors indicate more massive (luminous) stars (red = 
12-22 M©; yellow = 23-34 M Q ; green = 35-40; blue = 41-48 M ; cyan = 49-56 M Q ; white = ~56 M ). 
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Fig. 10. — The locations of the massive (>12 Mg) young stars in the vicinity of a star forming knot in our WFPC2 field north of the 
galaxy plotted on the SDSS color image (left) and our F606W WFPC2 image (right). Annuli mark our radial bins to search for a correlation 
between age and distance from knot center. Colors are the same as in Figure [9] 



